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Abstract 

The HARP and NA61/SHINE hadroproduction experiments as well as their implications for neutrino physics 
are discussed. HARP measurements have already been used for predictions of neutrino beams in K2K and Mini- 
BooNE/SciBooNE experiments and are also being used to improve the atmospheric neutrino flux predictions and to 
help in the optimization of neutrino factory and super-beam designs. First measurements released recently by the 
NA61/SHINE experiment are of significant importance for a precise prediction of the J-PARC neutrino beam used for 
the T2K experiment. Both HARP and NA61/SHINE experiments provide also a large amount of input for validation 
and tuning of hadron production models in Monte-Carlo generators. 

Keywords: hadron production measurements, accelerator neutrino beams, atmospheric neutrino fluxes, hadron 
production models 



1. The HARP experiment 

The HARP experiment 1 1 1 at the CERN PS was de- 
signed to make measurements of hadron yields from a 
large range of nuclear targets from hydrogen to lead and 
for incident particle momenta from 1.5 to 15 GeV/c. 
The main motivations are the measurement of pion 
yields for a quantitative design of the proton driver of a 
future neutrino factory and a super-beam O, a substan- 
tial improvement in the calculation of the atmospheric 
neutrino flux 1 3 1 and the measurement of particle yields 
as input for the flux calculation in accelerator neutrino 
experiments d, such as K2K fT.'Sl, MiniBooNE ||7| 
and SciBooNE (H. In addition to these specific aims, 
the data provided by HARP are valuable for validation 
and tuning of hadron production models used in simu- 
lation programs. The HARP experiment is described in 
more detail below in order to illustrate some common 
features of modern hadron production experiments. 

To provide a large angular and momentum coverage 
of the produced charged particles the HARP experiment 
makes use of a large-acceptance spectrometer consist- 
ing of forward and large-angle detection systems. A de- 
tailed description of the experimental apparatus can be 



found in Ref. 1 1 1. The forward spectrometer — based 
on large area drift chambers (9] and a dipole magnet 
complemented by a set of detectors for particle iden- 
tification (PID): a time-of-flight waU IM (TOFW), a 
large Cherenkov detector (CHE) and an electromag- 
netic calorimeter — covers polar angles up to 250 mrad 
which is well matched to the angular range of inter- 
est for the measurement of hadron production to calcu- 
late the properties of conventional neutrino beams. The 
large-angle spectrometer — based on a Time Projec- 
tion Chamber (TPC) located inside a solenoidal magnet 
— has a large acceptance in the momentum and angu- 
lar range for the pions relevant to the production of the 
muons in a neutrino factory. It covers the large majority 
of the pions accepted in the focusing system of a typical 
design. The neutrino beam of a neutrino factory orig- 
inates from the decay of muons which are in turn the 
decay products of pions. 

A full set of data collected by the HARP experiment 
with solid thin (5% of nuclear interaction length, Ai) and 
cryogenic targets have been analyzed and published (TTJ 
[l2l[L3, 14, 15l[l6l|T7l[l8l J9, 20, 21, 22, 23|. Those 
results cover all the physics subjects discussed above. 



Boris A. Popov / Nuclear Physics B Proceedings Supplement 00 (2012) 7-[S| 



2 




p [Gev/c] P [GeV/c] p [GeV/c] 



Figure 1: Measurements of the double-differential production cross-sections of n'^ (open circles) and n~ (closed circles) from 12 GeV/c protons 
on carbon (left), N2 (middle) and O2 (right) targets as a function of pion momentum, p, in bins of pion polar angle, 6, in the laboratory frame. 
Different panels show different angular intervals. The error bars shown include statistical errors and all (diagonal) systematic errors. The curves 
show the SW parametrization of Eq. |3j with parameters given in Ref. 1 14 1 (solid line for n'^ and dashed line for 



1.1. Results obtained with the HARP forward spectrom- 
eter 

A detailed description of established experimental 
techniques for the data analysis in the HARP forward 
spectrometer can be found in Refs. (TT] |T2| [131 lEl- 
Only a brief summary is given here. 

The absolutely normalized double-differential cross- 
section for a process like p + Target + X can be 
expressed in bins of pion kinematic variables in the lab- 
oratory frame (momentum, pj^, and polar angle, Oj^), as 

^^{p,,e,) = -A^ . -i— . ^.N^\p^^e,) , (1) 

dpdO. Np^pt ApAQ NpQi 

where is the cross-section in cm^/(GeV/c)/srad 
for each (pj^, 6j^) bin covered in the analysis; ^ is the 
reciprocal of the number density of target nuclei; t is 
the thickness of the target along the beam direction; Ap 
and AQ are the bin sizes in momentum and solid angle 

{Ap = Pniax - Pmin\ AQ = InicOSiOmin) " COSiOmax)))'. 

A^pot is the number of protons on target after event se- 
lection cuts; N^^(Pj^,Oj^) is the yield of positive pions 
in bins of true momentum and angle in the laboratory 
frame. Eq. ([T]) can be generalized to give the inclusive 
cross-section for a particle of type a 

where reconstructed quantities are marked with a 
prime and M~l , is the inverse of a matrix which 

^ puap 6 a 

fully describes the migrations between bins of true and 



reconstructed quantities, namely: lab frame momentum, 
p, lab frame angle, 0, and particle type, a. 

There is a background associated with beam protons 
interacting in materials other than the nuclear target 
(parts of the detector, air, etc.). These events are sub- 
tracted by using data collected without the nuclear tar- 
get in place after normalization to the same number of 
protons on target. This procedure is referred to as the 
'empty target subtraction'. 

The event selection is performed in the following 
way: a good event is required to have a single, well re- 
constructed and identified beam particle impinging on 
the nuclear target. A downstream trigger in the forward 
trigger plane (FTP) is also required to record the event, 
necessitating an additional set of unbiased, pre-scaled 
triggers for absolute normalization of the cross- section. 
These pre-scaled triggers (e.g. 1/64) are subject to ex- 
actly the same selection criteria for a 'good' beam par- 
ticle as the event triggers allowing the efficiencies of 
the selection to cancel, thus adding no additional sys- 
tematic uncertainty to the absolute normalization of the 
result. Secondary track selection criteria have been op- 
timized to ensure the quality of the momentum recon- 
struction as well as a clean time-of-flight measurement 
while maintaining high reconstruction and particle iden- 
tification efficiencies |[T2l[T3l . 

The first HARP physics publication |11| reported 
measurements of the tt^ production cross- section from 
a thin 5% Ai aluminum target at 12.9 GeV/c proton mo- 
mentum. This corresponds to the energy of the KEK 
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Figure 2: The dependence on the atomic number A of the forward proton production yields in p-A and 7r--A (A = Be, C, Al, Cu, Sn, Ta, Pb) 
interactions averaged over two forward angular regions (0.05 rad < 6 < 0.15 rad and 0.15 rad < 6 < 0.25 rad) and four momentum regions 
(0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 2.5 GeV/c, 2.5 GeV/c <p<3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for the four different 
incoming beam momenta (3, 5, 8 and 12 GeV/c). 



PS and the target material used by the K2K experiment. 
The results obtained in Ref. 1 11] were subsequently ap- 
plied to the final K2K neutrino oscillation analysis |6|, 
allowing a significant reduction of the dominant system- 
atic error associated with the calculation of the so-called 
far-to-near ratio from 5.1% to 2.9% (see 1 11 1 and | 6 | for 
a detailed discussion) and thus an increased K2K sensi- 
tivity to the oscillation signal. 

The next HARP goal was to contribute to the un- 
derstanding of the MiniBooNE and SciBooNE neutrino 
fluxes. They are both produced by the Booster Neutrino 
Beam at Fermilab which originates from protons accel- 
erated to 8.9 GeV/c by the booster before being collided 
against a beryllium target. A fundamental input for the 
calculation of the resulting flux is the measurement of 
the TT^ cross-sections from a thin 5% Ai beryllium target 
at exactly 8.9 GeV/c proton momentum. 

These double-difl'erential cross- sections were mea- 
sured in the kinematic range from 0.75 GeV/c < pj: < 
6.5 GeV/c and 0.030 rad < 6^ < 0.210 rad, subdivided 
into 13 momentum and 6 angular bins HTSl . A typi- 
cal total uncertainty of 9.8% on the double-difl'erential 
cross-section values and a 4.9% uncertainty on the total 
integrated cross-section are obtained. These HARP re- 
sults have been used for neutrino flux predictions [1241 
in the MiniBooNE |25 1 and SciBooNE experiments O. 

Sanford and Wang 1261 have developed an empiri- 
cal parametrization for describing the production cross- 
sections of mesons in proton-nucleus interactions (e.g. 
p-hA TT^ -\-X). This Sanford- Wang (SW) parametriza- 



tion has the functional form: 

^-^(P, e) = exp[ci-C3^-C6^(/7-C7/?b cos^^ e)]p'Hl- — ) , (3) 
dpdQ p^ p^ 

where X denotes any system of other particles in the 
final state; p\y is the proton beam momentum in GeV/c; 
p and are the n'^ momentum and angle in units of 
GeV/c and radians, respectively; d^cr/(dpdD) is ex- 
pressed in units of mb/(GeV/c sr); dfl = 2n d(cos6); 
and the parameters ci, . . . , cg are obtained from fits to 
meson production data. HARP data were fitted using 
this empirical SW parametrization. This is a useful tool 
to compare and combine diff'erent data sets. A global 
SW parametrization for forward production of charged 
pions as an approximation of all the studied datasets is 
provided in Ref. 1221 . It can serve as a tool for quick 
yields estimates. 

The next HARP analysis was devoted to the measure- 
ment of the double-diff'erential production cross-section 
of 71- in the collision of 12 GeV/c protons (see Fig.[T]) 
and pions with a thin 5% Ai carbon target 1 14|. These 
measurements are important for a precise calculation of 
the atmospheric neutrino flux f3 1 and for a prediction of 
the development of extended air showers. Simulations 
predict that collisions of protons with a carbon target 
are very similar to proton interactions with air. This hy- 
pothesis could be directly tested with the HARP data. 
Measurements with cryogenic O2 and N2 targets |15|, 
also shown in Fig. [T] confirm that p-C data can indeed 
be used to predict pion production in P-O2 and P-N2 in- 
teractions. 
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It is important to emphasize that a systematic cam- 
paign of measurements was performed by HARP: all 
thin target data taken with pion beams are analyzed and 
published now 1761 . Similar results were also obtained 
for incoming protons [ 22 1 . Moreover, results on forward 
proton production in p-A and 7i--A interactions were 
published recently 1231 . Fig.|2]shows the dependence on 
the atomic number A of the forward proton yields in p- 
A and n--A (A = Be, C, Al, Cu, Sn, Ta, Pb) interactions 
averaged over two forward angular regions (0.05 rad < 
e < 0.15 rad and 0.15 rad < 6> < 0.25 rad) and four 
momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 
1.5 GeV/c < p < 2.5 GeV/c, 2.5 GeV/c < p < 
3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for 
the four different incoming beam momenta (3, 5, 8 and 
12 GeV/c). 

The advantage of all these measurements is that they 
were performed with the same detector, thus related sys- 
tematic uncertainties are minimized. 

1.2. Results obtained with the HARP large-angle spec- 
trometer 
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Figure 3: The dependence on the beam momentum of the tt^ (left) 
and n~ (right) production yields in p-A (A = Be, C, Al, Cu, Sn, Ta, Pb) 
interactions averaged over the forward angular region (0.350 rad < 
< 1.550 rad) and momentum region 100 MeV/c < p < 700 MeV/c. 
The results are given in arbitrary units, with a consistent scale between 
the left and right panel. Data points for different target nuclei and 
equal momenta are slightly shifted horizontally with respect to each 
other to increase the visibility. 

First results on the measurements of the double- 
differential cross-section for the production of charged 
pions in p-A collisions emitted at large angles from the 
incoming beam direction were obtained for the full set 
of available nuclear targets using an initial part of the 
accelerator spill to avoid distortions in the TPC ITtHTsI 
[T9ll . The pions were produced by proton beams in a 
momentum range from 3 GeV/c to 12 GeV/c hitting a 
target with a thickness of 5% Xi. The angular and mo- 
mentum range covered by the experiment (100 MeV/c < 



p < 800 MeV/c and 0.35 rad < 6> < 2.15 rad) is of 
particular importance for the design of a neutrino fac- 
tory. Track recognition, momentum determination and 
particle identification were all performed based on the 
measurements made with the TPC (see {YT\ for more 
details). 

After the development of a special algorithm to cor- 
rect for the TPC dynamic distortions 1 27 1 and validation 
of the TPC performance with benchmarks based on real 
data |[28]| , full spill data were analyzed and published for 
incoming protons [W| and pions |21|. Fig. [3] illustrates 
the dependences of the measured charged pion yields on 
the beam momentum in p-A interactions averaged over 
the forward angular region (0.350 rad <6< 1.550 rad) 
and momentum region 100 MeV/c < p < 700 MeV/c. 

An additional comparison was also performed on the 
pion yields measured with the short (5% X\) and long 
(100% Al) nuclear targets |29 |. 

All HARP measurements described above were com- 
pared with predictions of Monte-Carlo models available 
within GEANT4 {M and MARS [31], as well as with 
the FLUKA (321 and GiBUU ES models. Some mod- 
els provide a reasonable description of HARP data in 
some regions (see e.g. [T4l[T5l[ra[T7l[T8l[T9ll2Ql[2n 
l23l and technical notes 1341 ), while there is no model 
which would describe all aspects of the data. 

Tables with HARP results in text format are available 
from the DURHAM database (H. 

2. The NA61/SHINE experiment 

The physics program of the NA61 or SHINE (where 
SHINE = SPS Heavy Ion and Neutrino Experiment) 
experiment at the CERN SPS 1361 consists of three 
main subjects. In the first stage of data taking measure- 
ments of hadron production in hadron-nucleus interac- 
tions needed for neutrino (T2K 1371 ) and cosmic-ray 
(Pierre Auger |38| and KASCADE f39|) experiments 
were performed. Later stages of the NA61 experiment 
study hadron production in proton-proton and proton- 
nucleus interactions with high statistics needed for a 
better understanding of high pt production at SPS en- 
ergies. Finally, the energy dependence of hadron pro- 
duction properties will be measured in p-p, p-Pb inter- 
actions as well as in nucleus-nucleus collisions, with the 
aim to identify the properties of the onset of deconfine- 
ment and find evidence for the critical point of strongly 
interacting matter. 

The NA61/SHINE apparatus is a large acceptance 
spectrometer at the CERN SPS for the study of the 
hadronic final states produced in interactions of vari- 
ous beam particles (tt, p, ions) with a variety of fixed 
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Figure 4: Examples of PID capabilities of the NA61/SHINE spectrometer for positively charged particles. (Left) Specific energy loss in the TPCs 
as a function of momentum. Curves show parameterizations of the mean dE/dx calculated for diff'erent particle species. (Middle) Mass squared, 
derived from the ToF-F measurement and the fitted path length and momentum, versus momentum. The lines show the expected mass squared 
values for diff'erent particles. (Right) Example of two-dimensional m^-dE/dx plot for the momentum range 2-3 GeV/c. Four clear accumulations 
corresponding to positrons, pions, kaons and protons are observed. The combination of both and dE/dx measurements provides close to 100% 
purity in the pion selection over the whole momentum range. 



targets at the SPS energies. The main components of 
the current detector were constructed and used by the 
NA49 experiment |40|. The main tracking devices are 
four large volume TPCs. Two of them, the vertex TPCs 
(VTPC-1 and VTPC-2), are located in the magnetic 
field of two super-conducting dipole magnets (maxi- 
mum bending power of 9 Tm) and two others (MTPC-L 
and MTPC-R) are positioned downstream of the mag- 
nets symmetrically with respect to the beam line. One 
additional small TPC, the so-called gap TPC (GTPC), 
is installed on the beam axis between the vertex TPCs. 
The setup is supplemented by time of flight detector 
arrays two of which (ToF-L/R) were inherited from 
NA49 and can provide a time measurement resolution 
of (Ttof < 90 ps. A new forward time of flight de- 
tector (ToF-F) with a time resolution of -120 ps was 
constructed in order to extend the acceptance of the 
NA61/SHINE setup for pion and kaon identification as 
required for the T2K measurements. The particle identi- 
fication in NA61 is based on the diff'erential energy loss 
dE/dx measured in the TPCs combined with the mass 
squared measurements based on the ToF information, 
see Fig.|4] 

Two carbon (isotropic graphite) targets were used 
for T2K-related measurements: 1) a 2 cm-long target 
(about 4% Ai) with density p = 1.84 g/cm^, the so-called 
thin target; 2) a 90 cm long cylinder of 2.6 cm diame- 
ter (about 1.9 Ai), the so-called T2K replica target with 
density p = 1.83 g/cm^. 

2.7. First thin-target results 

The NA61/SHINE experiment was approved at 
CERN in June 2007. The first pilot run was performed 
during October 2007. In total about 670 k events with 
the thin target, 230 k events with the T2K replica tar- 



get and 80 k events without target (empty target events) 
were registered during the 2007 run. Using these data 
interaction cross sections and charged pion spectra in 
p-C interactions at 31 GeV/c were first measured |41 1. 
Such measurements are required to improve predictions 
of the neutrino flux for the T2K long baseline neutrino 
oscillation experiment in Japan |[37]| . 

For normalization and cross section measurements 
we adopted the same procedure as the one developed 
by the NA49 Collaboration t43l . The measured in- 
elastic and production cross sections in p-C interac- 
tions at 31 GeV/c |44| are 257.2+1.9+8.9 mb and 
229.3+1.9+9.0 mb, respectively. 

Crucial for this analysis is the identification of the 
produced charged pions. Depending on the momentum 
interval, diff'erent approaches have been adopted, which 
lead also to diff'erent track selection criteria. The task 
is facilitated for the negatively charged pions, by the 
observation that more than 90% of primary negatively 
charged particles produced in p-C interactions at this en- 
ergy are tt", and thus the analysis of tt" spectra can also 
be carried out without additional particle identification. 

Three analysis methods were developed to obtain 
charged pion spectra. These are: 1) analysis of tt" 
mesons via measurements of negatively charged par- 
ticles (h~ analysis)', 2) analysis of tt^ and tt" mesons 
identified via dE/dx measurements in the TPCs (dE/dx 
analysis at low momentum) and 3) analysis of tt^ and 
TT" mesons identified via time-of-flight and dE/dx mea- 
surements in the ToF-F and TPCs, respectively (tof - 
dE/dx analysis). Each analysis yields fully corrected 
pion spectra with independently calculated statistical 
and systematic errors. The spectra were compared in 
overlapping phase- space domains to check their con- 
sistency. Complementary domains were combined to 
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Figure 5: Laboratory momentum distributions of ti^ mesons pro- 
duced in production p-C interactions at 31 GeV/c in different inter- 
vals of polar angle {0). The spectra are normalized to the mean tt^ 
multiplicity in all production p-C interactions. Error bars indicate sta- 
tistical and systematic uncertainties added in quadrature. The overall 
uncertainty (2.3%) due to the normalization procedure is not shown. 
Predictions of hadron production models, FLUKA2008 (solid line) 
and VENUS4.12 (dotted line) are also indicated. 



reach maximum acceptance. 

The agreement between spectra obtained by different 
methods is, in general, better than 10%. In order to ob- 
tain the final spectra consisting of statistically uncorre- 
cted points the measurement with the smallest total er- 
ror was selected. 

Inclusive production cross sections for negatively and 
positively charged pions are presented as a function 
of laboratory momentum in 10 intervals of the lab- 
oratory polar angle covering the range from up to 
420 mrad |41 1. The final spectra for tt^ are plotted in 
Fig.|5] For the purpose of a comparison of the data with 
model predictions the spectra were normalized to the 
mean 7i" multiplicity in all production interactions. This 
avoids uncertainties due to the diff'erent treatment of 
quasi-elastic interactions in models as well as problems 



due to the absence of predictions for inclusive cross sec- 
tions. A clear advantage of NA61/SHINE is that it fully 
covers the kinematic phase-space of interest for T2K. 

It is interesting to compare the tt^ spectra in p-C 
interactions at 31 GeV/c to the predictions of event 
generators of hadronic interactions. Models that have 
been frequently used for the interpretation of cosmic 
ray data, i.e. VENUS4.12 (H, FLUKA2008 |32| and 
URQMDl.3.1 149J were selected. They are part of 
the CORSIKA ISOll framework for the simulation of air 
showers and are typically used to generate hadron-air 
interactions at energies below 80 GeV. In order to ensure 
that all relevant settings of the generators are identical 
to the ones used in air shower simulations, p-C interac- 
tions at 3 1 GeV/c were simulated within CORSIKA in 
the so-called interaction test mode. The results are pre- 
sented in Fig. [5] More model comparisons can be found 
in the technical notes L51J . 




p [GeV/c] 



Figure 6: Comparison of measured spectra with model predic- 
tions. Distributions are normalized to the mean multiplicity in all 
production p-i-C interactions. The vertical error bars on the data points 
show the total (stat. and syst.) uncertainty. The horizontal error bars 
indicate the bin size in momentum. 



Recently, the tof - dEjdx analysis was applied to the 
same set of 2007 data in order to perform measurements 
of diff'erential production cross sections of positively 
charged kaons in p-^C interactions at 31 GeV/c |42|. 
The corresponding results are shown in Fig. |6] together 
with model predictions. The knowledge of kaon pro- 
duction is crucial for precisely predicting the intrinsic 
electron neutrino component and the high energy tail of 
the T2K beam. 

Moreover, preliminary results for protons were also 
obtained ||45l|46l, see Fig. [7] The analysis of neutral 
strange particle production using V^-like topology of 
their decays is currently ongoing. First preliminary re- 
sults can be found in WT\ . 



Boris A. Popov / Nuclear Physics B Proceedings Supplement 00 (2012) 7-[S| 



7 



0.01 



0.04 
0.02 



0<6<20 mrad 


20<9<40 mrad 


40<e<60 mrad 


: v\ 60<e<100 mrad ; 

i ' \ - 


- 100<e<140 mrad - 


ij^ 140<e< 180 mrad - 

1° V - 


180<e<240 mrad 

j^^^^^^^^^^^^^ NA61 Preliminary 


240<e<300 mrad 

* preliminary protons from 2007 

-FLUKA2011 

- VENUS4.12 



5 10 15 20 2 5 10 15 20 25 

p [GeV/c] p [GeV/c] 



Figure 7: Laboratory momentum distributions of protons produced 
in p-C interactions at 31 GeV/c in different intervals of polar an- 
gle {0). The spectra are normalized to the mean proton multiplic- 
ity in all production p-C interactions. Error bars indicate statistical 
and systematic uncertainties added in quadrature. The overall uncer- 
tainty (2.3%) due to the normalization procedure is not shown. Pre- 
dictions of hadron production models, FLUKA2011 (solid line) and 
VENUS4.12 (dashed line) are also indicated. 

2.2. First replica-target results 

First measurements of pion emission from the T2K 
replica target were also performed |52|. 

The main motivation for measurements of hadron 
emission from a replica of the T2K target is to reduce 
the systematic uncertainties on the prediction of the ini- 
tial neutrino flux originating from re-interactions in the 
long target. 

The long-target analysis presented here uses the low- 
statistics data collected in 2007 I 53]| . It however sets 
the ground for the ongoing analysis of high- statistics 
NA61/SHINE data with the repHca of the T2K target. 
It demonstrates that high-quality long-target data were 
successfully taken with the NA61/SHINE apparatus for 
T2K, and that such data can efl'ectively be used to con- 
strain the T2K neutrino flux predictions. A comparison 
of neutrino flux predictions based on thin-target hadron 
production measurements and long-target hadron emis- 
sion measurements was performed as an illustration of 
the complete procedure, see Fig. [8] A very good agree- 
ment between these two complementary approaches is 
observed. It is important to stress that this is the first 
complete example of application of long-target data for 
neutrino flux predictions. 

The data presented here and in Refs. ||4T] |42l [53l al- 
ready provide important information used for improved 



calculation of the T2K neutrino flux. Meanwhile, a 
much larger data set with both the thin (4% A\) and the 
T2K replica carbon targets was recorded in 2009 (about 
6 M triggers with the thin target and 2 M triggers with 
the replica target) and 2010 (about 10 M triggers with 
the replica target) and is presently being analyzed. This 
will lead to results of higher precision for charged pions 
and kaons, protons, and A. The new data will allow 
a further significant reduction of the uncertainties in the 
prediction of the neutrino flux in the T2K experiment. 

3. Conclusions 

The HARP experiment has already made important 
contributions to the cross-section measurements rele- 
vant for neutrino experiments. 

First measurements of charged pion and kaon produc- 
tion in p-C interactions at 3 1 GeV/c released recently by 
NA61/SHINE are of significant importance and were al- 
ready employed for precise predictions of the J-PARC 
neutrino beam used for the first stage of the T2K exper- 
iment. 

Both experiments provide also a large amount of in- 
put for validation and tuning of hadron production mod- 
els in Monte-Carlo generators. 

Thus, hadron production experiments have already 
contributed to the recent advances in neutrino physics. 
Clearly, hadroproduction studies are a must for future 
precision neutrino experiments. 
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